Electron-hole pair condensation in graphene bilayer 
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We consider the pairing of electrons and holes due to their Coulomb attraction in two parallel, 
independently gated graphene layers, separated by a barrier. At weak coupling, there exist the 
BCS-like pair-condensed state. Despite the fact that electrons and holes behave like massless Dirac 
fermions, the problem of BCS-like electron-hole pairing in graphene bilayer turns out to be rather 
similar to that in usual coupled semiconductor quantum wells. The distinctions are due to Berry 
phase of electronic wave functions and different screening properties. We estimate values of the 
gap in one-particle excitation spectrum for different interlayer distances and carrier concentrations. 
Influence of disorder is discussed. At large enough dielectric susceptibility of surrounding medium, 
the weak coupling regime holds even at arbitrarily small carrier concentrations. Localized electron- 
hole pairs are absent in graphene, thus the behavior of the system versus coupling strength is 
cardinally different from usual BCS-BEC crossover. 

PACS numbers: 81.05.Uw, 73.21.Ac, 74.78.-w 



Recent progress in experimental technology allowed a 
fabrication of graphene, one atomic layer separated from 
graphite crystal [1-4]. A lot of theoretical attention to 
graphene has emerged due to peculiar properties of its 
band structure, consisting in linear dispersion of elec- 
tron energy near two inequivalent points of the Brillouin 
zone [5-7]. An electron wave function close to these 
points is well described by the two-dimensional Dirac 
equation for massless particles [6, 7] with the Fermi ve- 
locity vf ~ 10 6 m/s f» c/300 playing the role of effective 
"speed of light" [2] . Several peculiar transport phenom- 
ena in graphene has been discovered experimentally, e.g., 
anomalous quantum Hall effect [3] and minimal conduc- 
tivity [2, 8]. Unique properties of graphene, such as un- 
usually high mobility of charge carriers [1] and a phase 
coherent transport [9, 10], allow to propose it as a base 
of future nanoelectronic devices [11-13]. 

In the present work we consider the formation of the 
condensate of spatially separated electron-hole pairs in 
bilayer graphene structure when interlayer tunneling is 
negligible. A condensation and superfluidity of spatially 
separated electron-hole pairs in usual coupled semicon- 
ductor quantum wells (CQW) due to their Coulomb 
attraction has been proposed theoretically in [14]. A 
nondissipative motion of resulting pairs leads to appear- 
ance of persistent electric currents, flowing in two layers 
in opposite directions (contrary to 3D case, where phase 
fixation leads to a formation of excitonic insulator state 
[15]). In the present Letter we consider another physical 
realization of a two-dimensional electron-hole system — 
graphene bilayer. Its schematic setup is shown on Fig.l. 

Two parallel graphene sheets are separated by a dielec- 
tric layer of thickness D, large enough to neglect tunnel- 
ing between them. By applying the gate voltage V g be- 
tween graphene sheet and a gate electrode, isolated from 
it by dielectric layer, one can adjust the charge carrier 



concentration, fixing the chemical potential fi at any de- 
sired level [1] (in addition to electrical one, chemical or 
electrochemical [16] doping of graphene is also possible). 
Chemical potentials in either of graphene layers may be 
adjusted independently. We consider the case of equal 
densities, when in the top layer the chemical potential 
is /i > 0, and charge carriers are electrons, whereas in 
the bottom layer the chemical potential is —fj, < 0, and 
charge carriers are holes. At weak coupling conditions, 
the system is unstable with respect to Bardeen-Cooper- 
Schrieffer (BCS) interlayer pairing of electrons and holes 
due to their Coulomb attraction. 

Consider the effective Hamiltonian of the system, re- 
sponsible for pairing of electrons from the top graphene 
layer and holes from the bottom layer. The influence 
of remaining part of the total Hamiltonian, correspond- 
ing to electron and hole interactions within individual 
graphene layers, manifests itself via screening of inter- 
layer Coulomb interaction. The effective Hamiltonian 
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FIG. 1: Schematic set-up of a system for realization of a 
pairing of spatially separated electron and holes in graphene 
bilayer. In the right: chemical potential positions in two 
graphene layers, adjusted by gate voltages V g and — Vg. 
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can be presented in the form: 



Hq = ££ k (a+a k + K b k) + i 



V(q) cos y cos ya+ +q 6+ _ q 6 k2 a kl) 

ki,k2,q 



(i) 



where at and 6k are destruction operators for Dirac elec- 
tron and hole quasiparticles with in-plane momentum 
k, £k = S«F|k| — n is a quasiparticle energy measured 
from the Fermi level, V(q) is the potential of screened 
electron-hole interaction, q is a momentum transmitted 
via Coulomb scattering, ipi and if2 are scattering angles 
for electron and hole, and S is the area of a bilayer. The 
graphene-specific factor cos(c/?i/2) cos(c/?2/2) results after 
folding over components of a spinor electron wave func- 
tion. The pairing is completely degenerate on spin and 
valley quantum numbers of electron and hole and thus 
we ignore these degrees of freedom. 

Within the framework of BCS approach for a pairing 
of electrons and holes with opposite momenta [17], the 
Hamiltonian {T]) transforms into the form: 

H = ^^{a+a k + b+b k ) + ^ x 



x£v(q) 



1 + COS if 



°lc+q fo -k-q 6 -kak, 



(2) 



where ip is angle between k and k + q, i.e. scattering 
angle, equal for electron and hole. Here the factor (1 + 
cos</?)/2 means an overlap of initial and final electron 
states. It originates from Berry phase of electronic wave 
functions in graphene and have no analogue in CQW [18]. 

The Hamiltonian ([2]) can be diagonalized by Bo- 
golyubov transformation: 

ak = Mk«k + v k f3\, 6_ k = -"k/3_k - v k a£. 
Introducing usual notations 
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we derive the self-consistent gap equation 



1 

~4^2 



dq^(q) 



1 + cos ip Ak-i 



2i?k+q ' 



(3) 



where we turned from summation to integration on q. 
Note that Eq. ^ differs from the analogous one for CQW 
in two aspects. The first one is the linear dependence of 
£k on k instead of a quadratic one for CQW; however, 
near the Fermi energy it can be linearized, hence this 
difference is not essential. The second aspect is the pres- 
ence of the overlap factor (1 + cos<^)/2, which suppresses 
backscattering in graphene. At weak coupling this factor 



turns into unity due to a predominance of small scatter- 
ing angles in the integral in ([3]) , but at stronger coupling 
it can weaken the pairing. 

The main contribution to the integral in ([3]) comes from 
the region near to the Fermi energy, where the dynami- 
cally screened interlayer electron-hole interaction V(q, ui) 
is attractive. In random phase approximation, 
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(4) 



where v q = 2ire 2 /eq is the bare Coulomb interaction, e is 
the dielectric constant of a surrounding medium, xi an d 
X2 are dynamic polarizabilities within the top and bot- 
tom graphene sheets. In the case of equal densities, due 
to the particle-hole symmetry, polarizabilities are equal 
in both graphene sheets: xi — X2 — X- The equation 

1 - 2t, qX (q,^) +t, qX 2 (q,^)(l „ e ~2 qD) = (5) 

describes two branches w±(q) of a plasmon dispersion 
in the system, corresponding to in-phase and antiphase 
plasma oscillations [14, 19, 20]. 

The weak coupling (or BCS) regime takes place when 
the region of pairing near to the Fermi energy is narrow 
with respect to the Fermi energy itself. In such a case, 
the radial integration on £ = £k+q in the gap equation Q 
can be decoupled from the integration on the polar angle 
ip. We restrict the integration on £ to the region limited 
by the cutoff energy fuv, where the dynamically screened 
interaction (j4]) is attractive. As for the integration on tp, 
we perform it with the statically screened electron-hole 
attractive potential V^q) = V(q, 0). 

Assume that the main part of the integral in (UJ comes 
due to integrands with small q, not larger than q by order 
of magnitude. Let define q as satisfying the equation 
V(q) = V(Q)/2. The explicit form of V(q) is determined 
by the static polarizability [20-22] 



x(q,o) = - 



2nh 2 v 2 



(0) 



where a is the Thomas-Fermi screening length in 
graphene. Coefficients g s = g v = 2 arise due to spin 
and valley degeneracy of electron states in graphene. 

Let introduce the Fermi momentum fco = fJ-/hv-p and 
the dimensionless parameter 

2ehvF 1 

= — « 0.23 x e, 



2r B 



where the dimensionless Wigner-Seitz radius r s measures 
the ratio of the characteristic Coulomb energy of quan- 
tum system to its characteristic kinetic energy [20]. In 
2D semiconductor system r s increases with decreasing 
carrier density, but in graphene it is determined only by 
a dielectric constant e of surrounding medium. 

There are three characteristic distances in the system, 
namely a, D and the mean separation between charge 
carriers within each graphene layer I ~ 1/ko- A behavior 
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of the system depends on a relation between a, D and I. 
As can be shown, a/l — a and D/l — k Q D. Therefore, a 
behavior of the system is governed by two dimensionless 
parameters, a and fco-D. 

The characteristic momentum q can easily be found 
from ((4]) and (J6j> in two limiting cases: a -C fco-D and 
a fco-D. In the first one, an effective momentum cut- 
off occurs due to the factor exp(— qD) in the numerator 
of (TJI, thus we get q « 2/D. In the second one, it is 
determined by the Thomas-Fermi screening, and we get 
q ~ 4k /a. Both results can be written as 



4fcn 



D 



(7) 



The cutoff energy hu> is determined by a characteris- 
tic frequency of the lower branch of plasma oscillations 
and can be estimated as h& = huj-(q). To the first order 
in the electron-electron interaction, the dynamic polariz- 
ability at q — > and u> > vpq is [20, 22] 



x(q,w) 



AttH 2 



(8) 



In the case a <C fco-D, from §5§ and (HJ) we find = 
v F q{k D/a) 1/2 and oj + (q) = v F (2k q/a) 1 / 2 . Thus the 
cutoff energy is Hu — 2fi/(koDa) 1 ^ 2 . In the case a 3> 
koD, the approximate expression |5]) is inapplicable, be- 
cause the lower branch of a plasmon dispersion, formally 
found with it, falls into a single-particle excitation con- 
tinuum uj < vpq. Actually, in this case oj-(q) = Vpq 
and uj+(q) = vp^koq/a) 1 / 2 , and the cutoff energy is 
hHi = 4/j,/a. 

Electron-hole Cooper pairs have a size of order of l/q 
in the in-plane direction. The weak coupling requires a 
large pair size relative to mean separation between near- 
est pairs, i.e. ql <C 1. Using ([7]), we conclude, that the 
weak coupling regime occurs when at least one of param- 
eters a or fco-D is large with respect to unity. 

It is not of necessity in our case for the pairing to be 
s-wave. Due to spatial separation, Pauli exclusion prin- 
ciple does not impose any conditions on relative angular 
momentum as well as spins and valleys of paired parti- 
cles. To seek an /-wave solution of ([3]) we assume that 
Ak = Aexp(iii/3k) at |£| < hu> and Ak = at |£| > hui. 
The approximate solution is 



i 47ra 1 
A = 2tux> exp \ — — | 



(9) 



Here the dimensionless /-wave harmonic of V(k) (cor- 
rected by the Berry phase factor) is 



Vi = 



e- koUx (l + cos<p)/2 
x + 4/a + 4(1 - e - 2k ° Dx ) /a 2 x 



llip d<p, (10) 



and x = q/ka = 2sin(</?/2) in the integral. 

The condensation will lead to the /-wave pairing with 
the largest A and thus the largest Vi- We have found 



numerically that at any values of a and fco-D the s-wave 
pairing dominates, though at very weak coupling a differ- 
ence between Vi for different I becomes negligible. So we 
consider only s-wave pairing hereafter. Deriving asymp- 
totic expressions for Vi and substituting them into ([9]) we 
find the gap for different relations between a, fco-D and 
unity. 

At a <C fcg-D we get 



4// 



(koDa) 1 / 2 



exp 



tirk Q D 1 



fco-D 



(11) 



This expression is suitable for values of a, both small and 
large with respect to unity. 

At q » feoD we distinguish cases of small and in- 
termediate inter layer distances D. In the first case, 
fco-D «l«a, the gap is 



8u f 27ra 

A = — exp <^ ; — 

a F \ In 1 + a/2 



In the second case, 1 <C fco-D <C a, the gap is 



8 M 

A = — exp 

a 



2na 



\n(a/4k D) - 7 



(12) 



(13) 



where 7 w 0.577 is the Euler constant. 

We have considered the weak coupling regime, but 
what can occur at stronger coupling? In CQW, at 
T = 0, on increase of a coupling strength there ex- 
ist a crossover from BCS-like state to Bose-Einstein 
condensation (BEC) in a dilute gas of localized, non- 
overlapping electron-hole pairs (or excitons for quasi- 
cquilibrium state created after laser pumping). At T ^ 0, 
both BCS-like state and gas of local pairs are in super- 
fluid state below the temperature of Kosterlitz-Thouless 
transition to the normal state [14]. In graphene, there are 
no bound solutions for the Dirac problem of single elec- 
tron in attractive potential due to the absence of a gap in 
the energy spectrum. Similarly, as can be shown, there 
are no localized electron hole-pairs in a bilayer. Therefore 
a behavior of graphene electron-hole bilayer on increase 
of the coupling strength will be cardinally different from 
BCS-BEC crossover in CQW. Strong coupling regime will 
be studied in subsequent publication. Note, that in a per- 
pendicular magnetic field the existence of localized mag- 
netoexcitons makes BCS-BEC crossover possible [23, 24]. 

Consider the weak coupling conditions more closely at 
a < 1 or a ~ 1. Such values of a can be realized 
with, e.g., commonly used SiC>2 substrate (e ~ 4). In 
this case the coupling strength is determined only by a 
value of fco-D. The Fermi momentum fco is proportional 
to /i, which can be tuned from zero to maximal values of 
« 0.3 eV in electrically doped graphene [1, 2]. The weak 
coupling regime (fco-D 3> 1) can be achieved with rea- 
sonable carrier concentrations at any interlayer distance 
D > 100 A. On the other side, by tending /i to zero, 
one can always achieve a strong coupling regime. Thus, 
the whole transition from weak to strong coupling can be 
realized experimentally by changing the gate voltage. 
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The case a ^> 1 takes place at large values of the dielec- 
tric constant e of surrounding medium (at least e > 5), 
which can be achieved with, e.g., HfC"2 (e ~ 25). In this 
case, the weak coupling regime sustains even at fj, — ► 0, 
i.e. at arbitrarily small carrier concentrations, and the 
gap tends to zero as A oc /i according to (fl2|) . This pro- 
vides a remarkable contrast with CWQ, where a strong 
coupling regime occurs inevitably at vanishingly small 
carrier concentrations. 

With the typical /i « 0.1 eV and minimal reasonable 
interlayer distance D = 50 A, at e — 7 (when weak cou- 
pling approximation is still reliable) the expression Q12p 
gives A = 4 x 10 -6 eV, which is equivalent to the temper- 
ature of 0.05 K. The maximal value of A can be achieved 
at strong coupling, when a ~ koD ~ 1. In this case 
A ~ /i, up to hundreds of Kelvins. 

Let now estimate the influence of disorder. In the 
case of conventional phonon-mediated superconductiv- 
ity, a presence of magnetic impurities acts destructively 
on BCS state because of a different scattering of elec- 
trons with opposite spins, which form Cooper pairs. In 
our case any impurity localized in either graphcne layer 
acts destructively on pair condensate, since it scatters 
only one pair constituent. The BCS-like state maintains 
if the mean free path A exceeds the coherence length 
I a = Hvp/A (analogously to Ref. [25]). From the ex- 
pression for diffusive Boltzmann conductivity [26] we de- 
rive A w (p,/ev-p)fi c , where /i c is the carrier mobility in 
graphene. In dirty graphene samples /i c « 1000cm 2 /V-s 
at room temperature [1] and the corresponding mean 
free path is A « 10 nm at // = 0.1 eV. On the con- 
trary, for clean graphene at the temperature of liquid 
helium /i c « 10 6 cm 2 /V ■ s, so A « 10 (im. Comparing 
these quantities with estimations of I a , we conclude that 
weak-coupling BCS-like state can be realized only in very 
clean graphene samples, whereas a strong-coupling state 
can survive at rather strong disorder. 

The onset of electron- hole pairing may be observed ex- 
perimentally via drag effect peculiarities. It has been 
shown for CQW, that an occurrence of pair condensa- 
tion leads to a sharp increase of drag resistivity [27-29]. 
Another possibility is an observation of both stationary 
and non-stationary Josephson-like effects (see, e.g., [30] 
and references therein). Also the condensation modifies 
an electromagnetic response of the system. In particular, 
an application of in-plane magnetic field leads to a for- 
mation of persistent dipolar supercurrent, which can be 
detected directly [14, 30, 31]. 

In CQW, the serious obstacle to an occurrence of a 
BCS-like electron-hole pairing is an anisotropy of the 
hole band [27, 32]. However in graphene the particle-hole 
symmetry provides almost perfect matching in shape be- 
tween electron and hole Fermi lines. In this connection 
it is interesting to discuss a possible effect of the trigonal 
warping in graphene on BCS-like pairing in our system. 
The trigonal warping breaks the isotropy of Dirac spec- 



trum in graphene and causes a deviation of the Fermi 
line from a perfect circle towards a triangle-like shape 
[18]. The warping is negligible at low carrier concentra- 
tions, but becomes considerable at large enough fi. The 
triangle-like deviation has opposite orientations in two 
graphene valleys and therefore breaks the valley symme- 
try, leading to the fixation of a condensate structure in 
valley space. However, in a case of small trigonal warp- 
ing a condensate with paired electrons and holes from 
different valleys will have slightly larger energy, than a 
condensate with electrons and holed paired from same 
valleys. In such a situation, a two-gap state with new 
collective excitation modes can be formed, similarly to a 
superconductor with overlapping bands [33] . 

Two spatially separated bilayer graphene sheets can 
also be a candidate for a realization of the BCS-like pair 
condensation. In a perpendicular electric field an elec- 
tron spectrum in bilayer graphene acquires a tunable gap 
and energy dispersion has a quite unusual "Mexican hat" 
shape [34, 35]. After electrical doping a nontrivial Fermi 
surface in the shape of a ring will be formed [36]. In 
principle, the BCS-like pairing between electron and hole 
Fermi-rings is possible, leading to a formation of two en- 
ergy gaps — inside and outside of the ring. 

In conclusion, we have analyzed a possibility of the 
BCS-like pairing between spatially separated electrons 
and holes in two parallel graphene layers with negligi- 
ble interlayer tunneling. At weak coupling the problem 
of pairing is rather similar to that in CQW, except for 
several graphene-specific differences, i.e. Berry phase of 
electronic wave functions and different screening prop- 
erties. We have derived asymptotic expressions for the 
gap in the excitation spectrum at various characteris- 
tics of the system and estimated its numerical value un- 
der reasonable conditions. Estimations of an influence 
of disorder have been carried out. An appearance of the 
BCS-like electron-hole condensate can be observed exper- 
imentally via Coulomb drag measurements, by studying 
a Josephson-like effect or by probing an electromagnetic 
response of the system. The trigonal warping in the elec- 
tron energy dispersion can lead to a formation of a two- 
gap state. 

There is no localized electron-hole pairs in graphene 
bilayer due to absence of a gap in the energy spectrum. 
Therefore a behavior of the system on increase of the 
coupling strength is cardinally different from BCS-BEC 
crossover in CQW. Both weak and strong coupling condi- 
tions can be achieved experimentally. If a dielectric con- 
stant of surrounding medium is large enough, then the 
weak coupling regime sustains at arbitrarily small carrier 
concentrations, in contrast to a situation in CQW. 
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